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ABSTRACT 

Most gamma-ray bursts (GRBs) observed by the Swift satellite show an early steep 
decay phase (SDP) in their X-ray lightcurve, which is usually a smooth continuation 
of the prompt gamma-ray emission, strongly suggesting that it is its tail. However, the 
mechanism behind it is still not clear. The most popular model for this SDP is High 
Latitude Emission (HLE) , in which after the prompt emission from a (quasi-) spherical 
shell stops photons from increasingly large angles relative to the line of sight still reach 
the observer, with a smaller Doppler factor. This results in a simple relation between 
the temporal and spectral indexes, a = 2 + (3 where F^, oc t~°'i/~^ . While HLE is 
expected in many models for the prompt GRB emission, such as the popular internal 
shocks model, there are models in which it is not expected, such as sporadic magnetic 
reconnection events. Therefore, testing whether the SDP is consistent with HLE can 
help distinguish between different prompt emission models. In order to adequately 
address this question in a careful quantitative manner we develop a realistic self- 
consistent model for the prompt emission and its HLE tail, which can be used for 
combined temporal and spectral fits to GRB data that would provide strict tests for 
the HLE model. We model the prompt emission as the sum of its individual pulses 
with their HLE tails, where each pulse arises from an ultra-relativistic uniform thin 
spherical shell that emits isotropically in its own rest frame over a finite range of radii. 
Analytic expressions for the observed fiux density are obtained for the internal shock 
case with a Band function emission spectrum. We find that the observed instantaneous 
spectrum is also a Band function. Our model naturally produces, at least qualitatively, 
the observed spectral softening and steepening of the fiux decay as the peak photon 
energy sweeps across the observed energy range. The observed fiux during the SDP 
is initially dominated by the tail of the last pulse, but the tails of one or more earlier 
pulses can become dominant later on. A simple criterion is given for the dominant 
pulse at late times. The relation a = 2 + P holds also as P and a change in time. 
Modeling several overlapping pulses as a single wider pulse would over-predict the 
emission tail. 

Key words: Gamma-rays: bursts - methods: analytical. 



1 INTRODUCTION 

Before the launch of the Swift satellite (Gehrels et al. 2004), Gamma Ray burst (GRB) X-ray afterglows were detected 
at least several hours after the burst (Soffitta et al. 2004 and references therein). They typically displayed a power law decay 
~ — t~^'^ around their detection time (De Pasquale et al. 2006). Swifts ability to rapidly and autonomously slew when the 
Burst Alert Telescope (BAT, observing in the energy range 15 — 350 keV; Barthelmy et al. 2005) instrument detects a GRB 
enables it to point its other instruments - the X-Ray Telescope (XRT, observing in the energy range 0.2 — 10 keV; Burrows et 
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al. 2005a) and UV/Optical Telescope (UVOT, observing at wavelengths 170 — 650 nm, i.e. from the optical to the near UV; 
Roming et al. 2005) - toward the GRB within tens of seconds from the GRB trigger time. The XRT thus filled the observational 
gap between the end of the prompt emission and the beginning of the pie-Swift afterglow observations several hours later. 
It revealed a complex behaviour usually consisting of three phases, followed by most GRBs, and referred to as a canonical 
light curve (Nousek et al. 2006), consisting of three distinct power-law segments where K, oc an initial (at t < tbreak,i, 
with 300 s < tbreak.i < 500 s) Very steep decay with time t (with a power-law index 3 < ai < 5; see also Bartherlmy et 
al. 2005; Tagliaferri et al. 2005); a subsequent (at tbreak.i < t < ibreak,2, with 10"^ s < tbreak,2 ^ 10* s) very shallow decay 
(0.2 ^ 02 ^ 0.8); and a final steepening (at t > <:break,2) to the familiar pve-Swift power-law behaviour (1 ^ ^ 1.5). In 
many cases there are X-ray flares superimposed on this underlying smooth component (typically during the first two phases, 
at t < ibreak,2; Burrows et al. 2005b; Falcone et al. 2006; Krimm et al. 2007), and in some cases there is a later (at tj > tbreak,2) 
further steepening due to a jet. 

The third decay phase (F,y (x t^"^) is the afterglow emission that was observed before Swift and is well explained by 
synchrotron radiation from the forward shock that is driven into the external medium as the GRB ejecta are decelerated, 
where the energy in the afterglow shock is constant in time (no significant energy gains or losses). The plateau (or shallow 
decay) phase can be explained either by pre- Swift models or by later models that have been developed especially for this 
purpose (Nousek et al. 2006; Panaitescu et al. 2006; Granot 2007). It could be due to energy injection, either by a tail of 
decreasing Lorentz factors at the end of the ejection phase (Rees &; Mfiszaros 1998; Sari &; Mfiszaros 2000; Ramirez-Ruiz, 
Merloni &; Rees 2001; Granot & Kumar 2006) or by a relativistic wind produced by a long lasting source activity (Rees & 
Meszaros 2000; McFadyen et al. 2001; Lee & Ramirez- Ruiz 2002; Dai 2004; Ramirez- Ruiz 2004), by an increasing efficiency of 
X-ray afterglow emission due to time dependence of the shock microphysics parameters (Granot, Konigl &; Piran 2006), by a 
viewing angle slightly outside the region of prominent afterglow emission (Eichler & Granot 2006), by a contribution from the 
reverse-shock (Genet, Daigne & Mochkovitch 2007) or by a two component jet model (Peng, Konigl & Granot 2005; Granot, 
Konigl & Piran 2006). 

The steep decay phase is observed in most bursts, and is in the great majority of cases a smooth continuation of the 
prompt emission, both temporally and spectrally (O'Brien et al. 2006). This strongly suggests that it is the tail of the prompt 
emission. Several explanations for this phase have been suggested in the context of previously existing models (Tagliaferri et 
al. 2005; Nousek et al. 2006), such as emission from the hot cocoon in the coUapsar model (Mfezaros & Rees 2001, Ramirez- 
Ruiz et al. 2002). The most popular model, by far, is High Latitude Emission (HLE) originally referred to as emission from 
a "naked" GRB (Kumar and Panaitescu 2000a). In this model the prompt GRB emission is from a (quasi-) spherical shell, 
and after it turns off at some radius then photons keep reaching the observer from increasingly larger angles relative to the 
line of sight, due to the the added path length caused by the curvature of the emitting region. Such late arriving photons 
experience a smaller Doppler factor. This leads to a simple relation between the temporal and spectral indexes, a = 2 + f3 
where F,^{t) oc t~°'v~^ , that holds at late times when t — to 3> At, where to and At are the start time and width of the pulse, 
respectively. The steep decay phase also shows a softening of the spectrum with time (see Zhang et al. 2007 and references 
therein). 

The consistency of the steep decay phase with HLE has been studied by several authors (Nousek et al. 2006; Liang et 
al. 2006; Butler & Kocevski 2007; Zhang et al. 2007; Qin 2008). However, some simplifying assumptions were usually made, 
which may afi^ect the comparison between this model and the observations. One such assumption is the choice of the reference 
time to for the steep decay, especially when the prompt emission consists of several pulses. Liang et al. (2006) find that when 
assuming the HLE relation a = 2 -\- fi and fitting for to its derived value is consistent with the onset of the last pulse of the 
prompt emission (or of the individual spike or fiare whose tail is being fit). Zhang et al. (2007) find that the HLE cannot 
explain the steep decays accompanied by a spectral softening, but can explain the cases with no observed spectral evolution. 
Barniol Duran and Kumar (2008) find that only 20% of their sample is consistent with HLE. Butler & Kocevski (2007) find 
that for a (physically motivated) time independent soft X-ray absorption (fixed Nh) the spectrum during the steep decay 
phase, is much better fit b}' an intrinsic Band function spectrum (Band et al. 1993) rather than b}' a power-law, and that the 
peak photon energy shifts to lower energies with time. Qin (2008) finds that such a behavior can, at least qualitatively, be 
produced for a delta function emission in radius with a Band function spectrum. It is therefore still a largely open question 
whether the temporal and spectral properties of the steep decay are consistent with HLE. Moreover, it appear that a physically 
motivated model for the prompt emission with realistic assumptions about the emission (e.g. over a finite range of radii with a 
Band function emissivity) is needed in order to address this question in a more quantitative and fully self consistent manner. 

The nature of the prompt GRB emission is what ultimately determines the properties of its tail. HLE is expected only 
in models where the prompt emission is from a quasi-spherical shell and turns off rather abruptly at some finite radius (or 
lab frame time). The best example for this type of model is internal shocks (Rees k, Mfezaros 1994; Sari k, Piran 1997) 
where variability in the Lorentz factor of the relativistic GRB outflow causes faster shells of ejecta to collide with slower sells 
resulting in shocks going into the shell over a finite range of radii (typically A.R ~ R). On the other hands, there are models 
in which HLE is not expected, such as in the case of isolated sporadic magnetic reconnection events within a Poynting fiux 
dominated outfiow (e.g. Lyutikov &; Blandford 2003) in which each spike in the GRB light curve is from a distinct small and 
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well localized region. Therefore, testing whether the steep decay phase is consistent with HLE would help distinguish between 
these two types (or classes) of prompt GRB models. This can be an important step toward identifying the basic underlying 
mechanism for the prompt emission, which is still one the the most striking open questions in GRB research more than four 
decades after the discovery of GRBs. 

In order to address this question, we develop a model for the prompt and its HLE tail that is physically motivated, 
realistic, and easy to use (fully analytic in its simplest version) in global joint fits (to all of the available data at all times and 
photon energies) of the prompt GRB and its SDP tail. Such global fits can provide a stringent and fully self-consistent test 
of HLE model for the SDP in GRBs. 

The prompt emission is modeled as the sum of a finite number of pulses. Each pulse corresponds to a spike in the 
prompt GRB light curve and has its own HLE tail. An individual spike is modeled as arising from a thin uniform spherical 
relativistic shell that emits isotropically in its own rest frame over a finite range of radii, while the observed fiux is calculated 
by integrating over the equal arrival time surface (Granot, Piran & Sari 1999; Granot 2005; Granot, Cohen- Tanugi & do 
Couto e Silva 2008) of photons to a distant observer. Our model is particularly suitable for internal shocks, which we focus on 
in this paper. For the emitted spectrum we consider the phenomenological Band function, which provides a good fit the the 
prompt emission spectrum of the vast majority of GRBs. We point out that our model can also be used for X-ray fiares, which 
appear to have temporal and spectral properties similar to the spikes of the prompt GRB emission. The main text provides 
the most useful results in an easy to use form, while the full derivations of these results are provided in appendixes in order 
to help understand their origin and make it easier to extend or generalize our model. We stress here that our main aim is not 
necessarily to uniquely determine all of the model parameters, which may be subject to various degeneracies and may prove 
hard when fitting to real data, but instead to test whether our model can provide a good fit to the data for any set of physical 
parameters. While such a good fit would still not prove that the HLE must be at work, it would definitely support HLE as a 
viable and arguably most plausible model. Our model for an individual pulse is described in § [21 and results for the fiux in 
the case for internal shocks with a Band function spectrum are given in § O The dependence of a single pulse on the model 
parameters is then investigated in § [H while § [5] discusses how to combine several pulses in order obtain to the total prompt 
emission and its tail. Both are intended to help the reader when using our model to fit data, which is one of the main aims of 
our paper. Our conclusions are discussed in §[6l This paper describes in detail our theoretical model and its main properties, 
and stresses some important caveats that one should keep in mind when using it to fit data in order to test the HLE model. 
In subsequent work we intend to confront it with Swift BAT+XRT data. 



2 DESCRIPTION OF THE MODEL 

2.1 The Basic Physical Model 

We consider an ultra-relativistic (F 3> 1) thin (of width <C R/T^) spherical expanding shell that emits over a range 
of radii Ro ^ R ^ Rf = Ro + AiZ. The emission turns on at radius Ro and turns off at radius Rf > Ro- The Lorentz 
factor of the emitting shell is assumed to scale as a power- law with radius, = rg(i?/7?o)"'" where To = r{Ro). The 
emission is assumed to be isotropic in the comoving frame of the shell, and uniform over the shell, i.e. the comoving spectral 
luminosity depends only on the radius of the shell, L'^i = L'^i (i?). As the main purpose of this work is to check the consistency 
of the tail of the prompt emission with HLE, we need to model the prompt emission. We therefore use for the emission 
spectrum the phenomenological Band function (Band et al. 1993) spectrum that provides a good fit to the observed prompt 
emission spectrum of the vast majority of GRBs. In the following we mainly consider emission over a finite range of radii, 
Ro ^ R ^ Rf — Ro + A_R. The comoving luminosity is then: 

^-'=^o( — ) -^l — )' S{x)^e 1 <; (1) 

'2 

are the high and low energy slopes of the spectrum. For &i > — 1 > &2 the Band function has a peak in the r^Fi, spectrum, at 
X = 1, and therefore since S{x) is normalized such that S{x) = xS{x) = 1 at a; = 1, it will not affect normalization of vFi, at 
its peak. The two functional forms used in the band function are matched at u'f, = Xhv'p. The peak luminosity L'^, evolve as a 
power-law with radius, L'^i = L'q{R/Ro)'^ where L'q = L'^, (i?o) is a normalization factor. 

Throughout the paper, primed quantities are quantities measured in the comoving frame (i.e. the local rest frame of the 
emitting shell), unprimed quantities are measured either in the source rest frame (the lab frame, i.e. the cosmological frame 
of the GRB; this includes F, R, 6 and t) or the observer frame (this refers to observed quantities, such as F,^, v and T). 



Ro) ' oy^,-. I ^^^^ 

where v'^(R) = v'o{R/RqY is the frequency where v' L'^,{R) peaks, with u'q = u'p{Ro); Xb — (fei — &2)/(l + &i), while 6i and k 
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2.2 Calculating the Observed Flux 

The observer is assumed to be located at a distance from the source that is much larger than the source size, so that 
the angle extended by the source as seen by the observer is very small and the observer effectively at "infinity". In order to 
calculate the fiux density Fi, that reaches the observer at an observed time T we integrate the luminosity L'^, over the Equal 
Arrival Time Surface (EATS; see Figure [T|), i.e. the locus of points from which photons that are emitted by the shell at a 
radius R, angle 9 relative to the line of sight, and a lab frame time t, reach the observer simultaneously at an observed time 
T (for full derivation see Appendix [A}. 



2.3 Expected parameters values for internal shocks 

The internal shocks model is the most popular model for the prompt GRB emission. Moreover, our model is very suitable 
for internal shocks. Therefore, we consider it in the following. Here we calculate the scalings of the various quantities with 
radius, that are expected for the internal shocks model. First, when different shells (i.e. parts of the outflow with different 
Lorentz factors) collide, they are expected to be in the coasting phase, corresponding to m = 0. Moreover, for the simplest 
case of uniform shells, the strength of the shocks going into the two shells, as characterized by the relative upstream to 
downstream Lorentz factor, Tud, is expected to be roughly constant with radius while the shock are crossing the shells. The 
electrons are expected to be fast cooling, i.e. cool signiflcantly on a timescale much shorter than the shell crossing time of the 
shock, and therefore most of the emission is expected to arise from a thin cooling layer behind the shock. Therefore our thin 
shell approximation is expected to be reasonably valid. Admittedly, we use one emitting thin shell, corresponding to a single 
shock front, while the shock going into the other shell is not explicitly modeled. One could always model such a second shock 
by adding another thin emitting shell that turns on and off at the same radii {Ro and Rf, respectively) but has a slightly 
smaller or larger Lorentz factor. This will not introduce a big difference in the overall result, so for the sake of simplicity we 
do not include this here. 

Now we turn to flnd the expected scaling of L'^, and v'p with radius, under the assumption that the observed soft 
gamma-ray range is dominated by synchrotron emission. For fast cooling, the peak frequency v'p of the uFi, spectrum is 
u'^ ^ {eB'"i'^)/{2-nme(?) where dNe/d-ye oc 7"'' for 7e > 7™ where 7m = (p - 2)/(p - l)(ee/Ce)(mp/me)(r„d — 1), while 
ee is the fraction of the internal energy behind the shock in the power law distribution of the relativistic electron, and 
is the fraction of all electrons taking part in this power energy distribution (and an electron-proton plasma is assumed for 
the composition of the outflow). As mentioned above, Vud is expected to be roughly constant during the shell crossing (for 
roughly uniform colliding shells), and therefore 7™ would also be approximately constant, so that v'p oc B' . The magnetic 
fleld is expected to be predominantly normal to the radial direction, so that B' ~ B/F cx B for m = 0. Moreover, B oc R~^ 
is expected both for a magnetic fleld convected from the central source, and for a fleld generated at the shock that hold 
some constant fraction (es) of the internal energy behind the shock. Therefore, one expects the peak frequency to evolve 
as u'p oc R^^. We have also assumed L'^, oc (R/Ro)". For synchrotron emission L^/ oc NeB' oc i?" as the number 
of emitting electron is proportional to the radius, A''e oc R. Since the cooling break frequency scales as v'^ oc R, we have 
K'^ ~ K',m!,A^'m/'^'c)~^^^ OC i?\ implying a = 1. 

More generally (without specifying the emission mechanism) for roughly uniform shells with constant Tud both the rate 
at which particles cross the shock and the average energy per particle are constant with radius, implying a constant rate of 
internal energy generation [dE'^^^/dt' (x R'^), and therefore for fast cooling this also applies for the total comoving luminosity, 
L' ~ v'pL'^i oc i?", and therefore d + a = Q. This is indeed satisfled for synchrotron emission for which d — —1 and a = 1, and 
holds more generally for other emission mechanisms in the fast cooling regime. 

For now on the values m — Q and d = — 1 derived in this part will be used throughout the paper. However, since the 
expressions do not become much simpler by specifying the value of a, we leave a in the simpler expressions, and use the value 
of a = 1 for flgures only. In particular, all the flgures showing lightcurves in this paper use these parameter values, as well as 
the mean BATSE values for the Band function spectral slopes: bi — —0.25 and 62 — —1.25 (Preece et al., 2000). 



2.4 Relevant Times and Timescales 

A photon emitted from the source (at the origin) when the shell is ejected from it (i.e. at a lab frame time toj when the 
shell radius is i? = 0) arrives at the observer at an observer time Tej which can be thought of as the observed ejection time 
of the shell. We deflne To the initial radial time by T = Tcj + To being the time at which the flrst photons emitted reach 
the observer (that is, photons emitted at a radius -Ro along the line of sight). We also deflne T/ the flnal angular time by 
T — Tc]+ Tf being the time at which the last photons that are emitted along the line of sight (from Rf and 9 = 0) reach the 
observer. 

For a constant Lorentz factor with radius (m = 0), as expected for internal shocks, the expressions for To and Tf are 
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simple: 

We also define two normalized times (and their corresponding values at T/) that will be used in the following: 

r = i + f=^^, f,^i + r, = g = i + ^ , (3) 

where T = 1 (or T = 0) corresponds to the onset of the spike - the very first photon that reaches the observer (emitted at 
Ro on the line of sight). The main motivation for defining these two times is that they correspond to the two most natural 
choices for the zero to, T = corresponding to the ejection time of the shell, and T = corresponding to the onset of the 
spike in the lightcurve. The choice of the zero time is important for the definition of the temporal index in S I4.1I where we 
explore these two choices in detail. Moreover, it is more convenient to use T in some expressions and T in others. 



L',, =L'oS[ — ] RoSiR - Ro 



3 RESULTS FOR INTERNAL SHOCKS WITH A BAND FUNCTION SPECTRUM 

3.1 Emission from a single radius 

Before to turn to the more generic case of emission from a range of radii, we first consider the limiting case of emission 
from a single radius Rq. The peak frequency is then Vp = v'q, and the luminosity is 

which after some algebra (see appendix [A] for details, and in particular section [A3|l we obtain the flux: 

P-iT ^ 1) = iltilioT-^S (-^' 

4:71(11 \'^0 

where and z are the luminosity distance and cosmological redshift of the source, Lq = 2roLo and fo = 2roi^o/(l + z). 
Denoting Fs = Lo(l + z)/{A-Kd\) and using the explicit expression for the Band function (eq. [T]), the observed flux reads 

r'i-2(i.//.o)''ie(i+'i)[i-^''/-«l f < xt,VQ/v , 

(6) 



(4) 



(5) 



b2bi-b2„l + b2 



3.2 Emission from a finite range of radii 

Integrating the luminosity (eq. ^) over the Equal Arrival Time Surface (for details of the calculation see appendix [Al 
and in particular its section [A4)l leads to the following expression for the flux: 



2 + a 

min (T,Tf) - 1 



F^{f I) = Fof-^ 
where Fo = (1 + z)Lq/[{2 + a)A-Kd\]. This can be explicited as: 



1^0 



(7) 



Fo 



XbVo/l' <f<ff, 

f > max 



(8) 



l]x 



bl-b2 l + i)2 



Note that the observed function has exactly the same shape as the local spectral emissivity - a pure Band function. This 
occurs only for m — Q and d = — 1. 

In terms of number of photons A'^ per unit photon energy E, area A and observed normalized time T (which is simply 
equal to F„/hE), this can be expressed as 



dN 



dEdAdT 
where 



(S,f ^ 1) =f- 



a + 2 

min ( T, Tf ) - 1 



e ^ z ^ hi ~ 1)2 

(fei -62)''i"'"e-<''i-''^' 61-62 



(9) 



(10) 
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is the familiar Band function with a normalization constant -Bnorm, where z = {E/Eo)T = (1 + bi)x, while E = hv and 
E(, — hvo are the corresponding photon energies (the more common notation is asand — hi — 1 and /Saand = 62 — 1). 



4 PROPERTIES OF THE SINGLE PULSE EMISSION 

Now that we have derived the observed flux for a single emission episode (or single pulse in the light curve), we study 
its temporal and spectral behaviour for any radial width Ai? ^ of the emitting region. We remind the reader that we 
consider only internal shocks, and use the corresponding model parameter values (a = 1, m = and d = —1) for fast cooling 
synchrotron emission, with a Band function emission (and observed) spectrum (except in some cases where the discussion can 
stay general without much complication). Some of the results may not hold for more general parameter values of m or d, and 
we point this out when relevant. For all flgures showing lightcurves (throughout the whole paper), the panels or figures with 
a linear scale show -F^/Fmax where Fmax = Fv{Tf), while panels or figures with a logarithmic scale show F^/Fo where we 
remind the reader that Fo{a = !) = (! + z)Lo/{12ndl). All fi gures showing temporal evolution of parameters or lightcurves 
with a logarithmic time axis in this section use T, as this shows the early behaviour much more clearly than for T. 

From eq. for reasonable values of the parameters T/, 61, 62, u/vq, and a, the pulse peaks at T = Tej + T/ (T = T/). 
While this is generally the case, for some combinations of parameters (often involving relatively large values of T/) the pulse 
has a round peak and starts decaying before Tf. 

For T < 1, the Equal Arrival Time Surface (EATS) does not intersect the emission region and no photons reach the 
observer (its outermost radius Rl is smaller than Ro): F^{f < 1) = 0. When 1 f ^ T/ (iio ^ < Rf), the EATS 
intersects the emission region but does not yet encounter its outer edge (in particular the observed fiux is independent of the 
radial extension AR of the emission region); the fraction of the EATS within the emission region increases with time, as does 
the maximal angle Smax relative to the line of sight from which photons reach the observer, (^maxFo)^ = (T— 1). When T > Tf 
{Rl > Rf), the front part of the EATS is outside the emission region, and its parts inside the emission region are at increasing 
angles from the line of sight. In particular, photons reach the observers from Bmin ^ ^ ^ 6'max where (SminFo)^ — {T — Tf)TJ^ . 
Note that at T ^ T/, well into the tail of the pulse, Smax/^min ~ T/^^, so that the emission comes from a rather narrow 
range of angles 9, whose typical value increases linearly with T. Moreover, for T > Tf, the fiux ratio for two identical sets of 
emission parameters that differ only in their Tf (denoted by subscripts 1 and 2), is constant in time and equal to 

F,(f'>f/,i > 1) ff+^-l' F4f>ff,i^l) 2 + a ' 

The first ratio approaches AR2/AR1 for Ai?i,2 ^ -Ro, since this corresponds to the thin shell limit, while the overall emitted 
energy is proportional to A_R, since L'^,{R) ^ L'^,{Ro) is almost independent of R within the very thin emission region. For 
the second ratio, the denominator is the fiux for a delta function emission with radius, for which the total emitted energy 
is held fixed, and therefore the ratio approaches AR2/R0 ^ 1 in the limit of a thin emission region. The fact that the fiux 
ratio is constant in time at T > Tf holds only for m = and d = — 1, and means that the fiux at these late times (typically 
after the peak of the spike, which is usually at Tf) has the same time dependence regardless the width of the emitting region 
{AR). This can simplify the calculation of the fiux for a family of pulses that differ only in AR: one can calculate the fiux 
for AR = {Tf = 1) and apply it to f ^ Tf, multiplied by a factor [ff+'' - l]/(2 + a) for any value Ai? > {Tf > 1). 
Moreover, it is also sufficient to calculate the fiux for AR 00 and apply it to T < Tf (this holds much more generally; 
Granot, Cohen- Tanugi & do Couto e Silva 2008). 

Figure [2] shows light curves for a single pulse in both linear and logarithmic scales, for different values of the normalized 
frequency u/uq. The peak time is at Tf = 2 (equivalent to Tf — 1). The light curves sample the two parts of the Band 
function both before and after the peak time. The differences between the light curves for different frequencies refiect the 
spectral evolution, and in particular the evolution of the spectral break frequency Up. At higher observed frequencies u the 
change in the spectral and temporal indexes associated with the passages of Up occurs earlier. The shape of a pulse (left panel 
of figure [2]) can vary from being very spiky {dotted line) to a rounder peak {dot-dashed line), depending on the frequency. It 
may thus provide some latitude in the fitting of actual observed pulses. 

Figure [3] shows the dependence of the same pulse on Tf for three values of the normalized frequency u/fo (0.01, 0.1 and 
1). It is evident from the logarithmic scale figures that at T < T/ the fiux is independent of AR (and therefore of Tf), and 
that at T ^ Tf all the light curves have the same time dependence (i.e. their fiux ratio is constant in time). At any given time 
the spectrum is independent of AT? (this is valid only for m = and d = —1). The bottom right panel of this figure shows 
linear scale to help visualise a case where the peak of the pulse is before Tf. 

Figure m shows the dependence of the same pulse on the parameter a for three values of the normalized frequency u/vo 
(0.01, 0.1 and 1). We can see that, compared to the case for a = 1, when a increases the peak is at Tf and becomes sharper. 
When a decreases the pulse becomes larger, the slope for T ^ Tf becoming closer to zero up to a point where is is zero. For 
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values of a even smaller, the peak of the pulse occurs before Tf and becomes rounder; in this case at Tf only a sharp break 
is observed. 

4.1 Local temporal and spectral indexes 

It is natural to define the local values of the spectral and temporal indexes as the logarithmic derivatives of the flux density 
with respect to frequency and time, respectively. For the spectral index, there is no ambiguity and /3 = —d\ogFv/d\ogv. For 
the temporal index, however, we must choose a reference time, and the choice is not obvious, for this reason we consider two 
alternative deflnitions: Qej = — dlog _Fi,/dlog T, that uses the ejection time as the reference time, and aon = — dlog F^/dlog T 
that uses the onset of the spike as the reference time. The flgures in this subsection use the observed frequency u instead of its 
normalized value u/vo, in order to provide a more realistic example that could be at least qualitatively compared with data, 
and include the BAT and XRT energy ranges. For these figures we consider Eq = 2ro£'o/(l + z) = 300 keV, which could for 
example correspond to i5o = 1 keV, To = 300 and z = 1. 

Figure [5] shows the evolution of the temporal indexes Uej and aon during a pulse (See appendix [B] for the detailed 
evolution of the temporal and spectral slopes). The temporal index aej starts at very negative values and gradually increases, 
until at Tf it makes an abrupt jump to its value during the decaying part of the pulse, which is exactly 2 + /3 (see eq. |B7) 1. 
The temporal index Qon starts at early times, T ^ I, either at —1 for 7/ > and T < Tf, or from for Tf — » 0. Moreover, 
for < 7/ <C 1, don ~ —1 for T <Tf and Qon ~ for f/ < T <C 1 (see eqs. [B2| and |B4) 1. Note that when aon jumps from 
its negative value to a positive value at T = Tf (i.e. at the peak of the spike), it reaches the same function of T, independent 
of the time of the jump, Tf, and therefore the same function also holds for Tf = AR/Ro = (see eqs. |B2| and |B4| 1. At late 
times, T S> 1 and T > Tf, the HLE relation is approached, aon « 2 — 62. 

The left panel of figure [6] shows the evolution of 2 + /9 (where /3 is the spectral index) with the temporal indexes a^j and 
aon- The spectral index naturally softens (/3 increases with time), similar to what is typically observed (at least qualitatively), 
until it reaches — &2 at late times (T ^ XbVo/i' — 1). The change in /3 occurs earlier at higher photon energies. At T ^ Tf, 
aej =2 + 13 while aon only approaches 2 + /? at late times. 

In order to get a better idea of how the observed spectral index /3 is expected to behave in Swift XRT observations, we 
calculate its average values over the XRT energy range (0.2-10 keV). We define two average values, by integrating over either 
the frequency v or its logarithm log v: 



The middle panel of figure [6] shows the evolution of these two averages as well as the values of /3 at z^min = 0.2 keV, 
i^max/2 = 5 keV and //max = 10 keV. As expected, (/3)^ gives a larger weight to higher frequencies compared to (/9)iogi/, and 
its value it is usually very close to the spectral slope at !/max/2 (5 keV), except when the break frequency Vp of the Band 
spectrum passes through the XRT range, and the change in /3 within this range is the largest. Therefore, {/3)iog^ appears to 
better refiect the spectral slope measured over a finite frequency range. 



4.2 Spectrum 

The local spectral emissivity in the comoving frame is taken to be a Band function. We have seen previously that for 
the parameter values relevant for internal shocks (m = 0, d = — 1), the observed spectrum is also a pure Band function. 
This is evident in the right panel of figure [H which shows the temporal evolution of the observed spectrum in our model. It 
results from the fact that for these parameter values the observed peak frequency Vp is constant along the EATS. We have 
Vp/vo = Ep(T)/Eo = 1/r = 1/(1 + T) (see eq. l|A18p i which is independent of Tf. This behaviour is evident in the right 
panel of figure[6l where Ep/Eo is 1 at the onset of the spike (T = 0), Ep/Eo = 1/2 at the peak of the spike (T = Tf = 1), 
and Ep/Eo decreases roughly linearly with T at later times, during the tail of the pulse. 

5 COMBINING PULSES TO OBTAIN THE PROMPT EMISSION 

There is good observational evidence that the steep decay phase is the tail of the prompt emission (O'Brien et al. 2006). 
Within our model, the prompt emission is the sum over a finite number of pulses, and therefore the steep decay phase is the 
sum of their tails. In this section we provide examples of combining several pulses to model the prompt emission, and study 
the effect of varying the different pulse parameters. To this end, we start with a simple prompt emission model consisting 
of six pulses that are identical except for their ejection time Tcj (see Fig. [TK). Each pulse corresponds to a single emission 
episode of a particular shell that was ejected at Tcj^i (for i'th pulse), has an initial radial time To.i, and a final angular time of 
Tf^i. Then, we study the effect of changing the other model parameters one by one among the pulses. All lightcurves in this 





(12) 
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section are drawn against T, as the ejection time is different for each pulse (and then the definition of a T T would differ for 
each pulse). In Fig. [TJj the peak fiux fpcak is varied. Next, we vary To and/or Tf. In Fig. [3;, To is varied while Tf /To remains 
constant. In Fig. [Til, Tf and AR/ Ro vary while To and -Ro remain constant. In Fig. [7^, To and AR/ Rq vary while Tf and Rf 
remain constant. Each of these panels show the light curve in logarithmic scales, T = is set to the onset time of the first 
pulse, which means that Tcj,! = —To,!, thus showing the modeled prompt from a time close to what would be the trigger time 
for an observed burst. The red solid line represents the total prompt emission (the sum of all the pulses), while the black non 
solid lines are the individual underlying pulses. All the examples shown here of the prompt emission are for u/fo = 0.1. 

In the case of six equal pulses (Fig. [TK), later pulses appear to decay much more steeply just after their peak in a 
logarithmic scale with the zero time near the beginning of the first pulse. At very late times the relative contribution from 
the different pulses becomes almost the same. As the only parameter that varies between pulses is the ejection time, Toj, this 
change in slope must depend only on it. Noting that the temporal slope is q = — dlog Fi,/dlog T — aej/{i —Taj/T), and that 
the value of a^j just after the peak is independent of Tej (it depends only on Tf, see eq. |B6) 1. we can see that the value of a 
just after the peak scales as apeak = Q!ej,peak(l + T^j/Tf), since T = Tej + Tf is the time of the peak of the pulse. Since the 
pulses are equal they have the same Tf and a<;j,pcak, it is clear that Qpcak increases with Toj. At late times when T 2> Tcj, a 
approaches aej = 2 + 13. 

When varying T'peak while fixing the other parameters (see Fig. [Th), the relative fiux from each pulse at very late times 
is proportional to its Fpeak, so that the largest contribution is from the pulse with the largest Tpcak. 

At late times the observed fiux density of a single spike scales as T^ oc t'^"^ (see, e.g., eq. [8]). If at the peak time of 
the spike, which for simplicity is assumed here to be at T = Tcj + Tf (as is usually the case), the observed photon energy is 
at the high-energy part of the Band function, £ ^ = Eo{To/Tf)xi, or u ^ = i'o{To/Tf)xt, then (using eq. [8]) the fiux 
from the peak onwards is simply given by 

T.^.. (T > Tej + Tf) = F,.,peak ^^-^T^) ^ ' (13) 

while for E/E^ = v/v^ < 1 the expression is slightly more complicated, 

-^Y~\-i^+''l)(-/-o)i,T-T^,-Tf)/To y^, ^ r-Tcj ^ ToXbUo/v , 



^;y,peak 



(14) 

but the qualitative behaviour is still rather similar. Therefore, the fiux ratio of two pulses with ejection times Tcj^i < Tcj, 2 
and a comparable Ep{Tf) = {To/Tf)Eo (as is usually the case for different pulses in the prompt emission of the same GRB), 
at late times (T > max(Tcj,i + T/,i, Tcj,2 + T/,2) and T — Tcj, 2 3> Tcj, 2 — Tcj,i) is approximately 

^ (IllV^' , for min[T,.i,r,,2]>l and T - Tcj,2 » Tcj,2 - Tcj,i , (15) 

i'i/,2(i ) -t^pcak,2 \lf,2 J 

where /3 = — ^2 for v ^ while (3 is generally intermediate between — &2 and — bi for v <v,. 

Fig. [2t demonstrates this nicely for a series of six pulses with the same Tpcak but decreasing Tf , so that the later pulses 
with a smaller T/ decay faster and become sub-dominant at late times. At the latest times the first spike, which has the largest 
Tf, dominates the observed fiux in the tail emission. A similar behaviour is also seen in Fig. [Tli. In Fig. [7^ both Fpeak and Tf 
are the same between the different pulses, and therefore their tail fiuxes at late times are similar. In Fig. [Jh, To and T/ are 
varied while Tf /To is constant, and it can be seen that this corresponds to a rescaling of the pulse width (its typical duration) 
without effecting its shape. In Fig. [T]!, T/ and Tf /To are varied while To is constant, and this nicely demonstrates how the 
shape of the pulse depends on Tf /To- Typically, the rise time of a pulse is T/ — To while its decay time is Tf, so that the ratio 
of the rise and decay time is 1 — To/Tf. In Fig. [7^, To and T//T0 are varied while T/ is constant. In this case the rise time 
varies considerably between the different pulses while the decay timescale and the late time tail of the pulses are practically 
the same. This arises since the tail is dominated by emission from R ^ Rf, that in this case is very similar for all the pulses. 
Moreover, for the particular choice of parameters in Fig. [7^, where Ep{Tf) = {To/Tf)Eo and i?, = xtEp{Tf) remain constant 
for all the pulses, their late time tails have the same fiux normalization. This can be understood from eq. I|14p . where the fiux 
for f > ma.x{ff,XbVo/v) can be written as T^/F„,peak = {E / E^f'^-''^ exp[(6i - b2){E - E,) / E,]f*'^-'^ . 

Fig. [7f shows a more realistic example of the prompt emission, in which a larger number of model parameters is varied 
between the different pulses. This example contains only three pulses in order to be clearer. It can be seen that the fiux 
during the decaying phase is initially dominated by the last pulse just after its peak (T > 27 s), but the second peak becomes 
dominant (even if only by a small factor) as early as T ~ 37 s, and finally at T ~ 140 s the first pulse becomes the dominant 
one. This demonstrates that different pulses can dominate the observed fiux during the course of the steep decay phase. Which 
pulses would contribute more to the steep decay phase can be estimated according to their typical width (or duration), peak 
fiux, and peak time. The peak time is most important at the beginning of the steep decay phase, where the last spike always 
dominates just after its peak if its peak is above the fiux from the other spikes. Later on the relative contribution of the 
different spikes can be estimated according to eq. I|15p . Since the late time fiux scales as TpcakT^^*^ and usually < /3 < 2, 
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the power of Tf (which corresponds to the typical width of the spike) is higher than that of -fpcak, so that wider spikes tend 
to dominate over narrower spikes, even if their peak flux is somewhat lower. 

One should be very careful when fltting actual data with such a model. Fig. [8] shows what can happen if because of noisy 
data or coarse time bins, a prompt emission (red solid line) which is actually composed by several pulses (three, six or twelve 
in the cases shown; black non-solid lines) is fitted by a single broad pulse (green solid line). In this case the tail of the prompt 
emission can be significantly overestimated at late times, by a factor that tends to increase with the true number of underlying 
pulses. This can be understood by the simple example of comparing a single spike with A'' identical spikes with the same peak 
flux but a duration smaller by a factor of A'', for which the sum of their late time tail flux would be smaller than that of the 
single pulse by a factor of ~ A''^^''. However, in more realistic examples, the late time flux would often be dominated by the 
widest underlying pulse, so that its width would be more important than the total number of narrower underlying spikes. It 
is important to keep this effect in mind when confronting such a model with actual data. 



6 DISCUSSION AND CONCLUSIONS 

We have presented and explored a model for the prompt GRB emission and its high latitude emission (HLE) tail. This 
model is physically motivated and realistic: it consists of a finite number of emission episodes, each of which corresponds to 
a single spike in the prompt light curve, and is modeled by a relativistically expanding thin spherical uniform shell emitting 
isotropically in its own rest frame within a finite range of radii. Our model thus describes the prompt emission and the steep 
decay pahse as a whole from its very start to its late tail. Yet this model is easy to use (fully analytic in its simplest form 
described here), making it particularly suitable for detailed combined temporal and spectral global fits to the prompt GRB 
emission and the following steep decay phase (SDP). Such fits can provide a stricter test of the HLE model for the SDP 
compared to most previous models, since we use a single self-consistent model to fit both the prompt emission and the SDP, 
while most previous models fit only the SDP and are largely decoupled from the details of the prompt emission. Moreover, 
our model is also physically motivated, and more realistic than previous models. We have derived analytic expressions for the 
flux in the realistic case of a Band function spectrum (eqs. [7] and [8]), which consists of two power laws that smoothly join 
at some typical photon energy. 

The temporal evolution of the instantaneous values of the temporal {a) and spectral (/3) indexes for a single emission 
episode was studied, corresponding to a single observed pulse in the light curve. The definition of a is not unique as it depends 
on the choice of reference time. We explored two options for the reference time, either the ejection time (aej) or the onset 
time of the spike (aon), and found that for the former the HLE relation [a^j — 2 + (3) \s satisfied from immediately after the 
peak of the spike (T > Tj), while for the former it is only approached at late times (ofon ~ 2 + /3 at for T > Tf and T S> 1). 

We have intentionally chosen a simple model to describe the pulses, in order to reduce the number of free parameters. For 
a single emission episode (or pulse), in the most generic case there are ten free parameters: the power m = — 2dlog F/dlog i?, 
d = dlogv'p/dlogR, a = dlogL^, /dlogR, the normalization factor Fo (or Lo), three additional parameters for the Band 
function (the two spectral slopes, bi and 62, as well as the peak energy at the onset of the pulse Eo), the two timescales To 
and Tf, and the ejection time Toj. We have the general constraint AR > 0, which implies Tf = 1 + AR/Ro > 1. Focusing 
on the internal shocks model fixes some of these parameters: as the outflow is typically in the coasting phase, m = 0, while 
for synchrotron emission from fast cooling electrons d = —1 and a = 1. Since we expect AR/Ro = T/ ~ 1 we can flx T/ ~ 1 
(although a wider range, such as 0.2AR/Ro < 5, may still be considered as plausible). Fixing m, d, a, and Tf/To in this 
manner would leave only six free parameters. For a prompt emission with several pulses, one may be able in some cases to 
neglect the spectral evolution and use the same values of 61, 62, and -Eo for all the different pulses (or at least two of them, 
e.g. fei and 62), which leads to a total number of free parameter of 3(Afpuiscs + 1) (or 4A'puiscs 4- 2 if i5o cannot be flxed for all 
the pulses) for a burst with A^'puiscs pulses. 

The shape of the pulses in our model can vary considerably, from very spiky peaks to rounder ones, from a very sharp 
rise to shallower rise, and so on (see Figs.EJ-H). This can help reproduce some of the observed diversity in the shape of spike 
in the prompt light curve. This appears to be a promising feature of our model. However, we have an abrupt change in the 
temporal index at Tf , that usually corresponds to a sharp peak of the spike. This is caused by our model assumption that 
the emission abruptly shuts off at the outer emission radius Rf. Therefore, we also consider an alternative and more realistic 
assumption, which leads to a rounder peak for the spikes, where the emission more gradually turns off at i? > J?/. This is 
done by introducing and exponential turn-off with radius of the comoving spectral luminosity, L^/ {R) , and is examined in 
Appendix [O The more gradual the turn-off of the emission with radius the rounder the peak of the pulse in the light curve. 
This can help fit the observed variety of pulse shapes even better (at the cost of adding an additional free parameter). 

In the particular case of synchrotron emission from internal shocks, we find that the observed spectrum has the same 
shape as the emitted one, which in our case is modeled as a Band function. The observed peak photon energy of the Band 
function decreases with time, Ep(T) = Eo/T, naturally leading to a softening of the spectrum with time, similar to what is 
observed by Swift. Thus, our model can at least qualitatively reproduce the main temporal and spectral features observed by 
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Swift. The spectral index /3 evolves from its value below Ep {/3 = —bi) to its value above Ep (/3 = — &2), where the transition 
that corresponds to the passage of Ep through the observed energy band occurs at earlier times for higher observed photon 
energies (or frequencies). 

When modeling the prompt emission by combining several pulses, the SDP is initially dominated by the last pulse (just 
after its peak, if it is above the flux fro the other pulses), but can later be dominated by the tail of other pulses. The relative 
contribution of a pulse to the late time flux scales as ~ FpcakT^^^ . and therefore wider pulses (with a larger T/), and to a 
lesser extent pulses with a larger peak flux (Fpeak), tend to dominate the late time flux, deep into the SDP. Moreover, often 
the contribution to the total flux from the tails of several pulses can be comparable, so it cannot be adequately modeled 
using a single pulse model. Therefore, we caution here that modeling the steep decay phase using the HLE of a single pulse, 
F^ oc {T — Tre()~^^'^'^\ may lead to wrong conclusions, and all the more so if the reference time Tref is arbitrarily set to the 
GRB trigger time. Even if Tref is set to the onset time of the last spike, this may still be a bad approximation in many cases 
since (i) we find that aej = 2 + f3 (with Tref = Tcj) rather than tton (with Tief = Tcj + Tii, corresponding to the onset of the 
spike) while Qon approaches aej =2 + 13 only at late times well after the peak of the last pulse, and (ii) at such late times 
the flux often becomes dominated by the tails of earlier pulses. 

Our model can produce different shapes for the tail of the prompt emission, from close to a power law (which can have 
a different temporal index than its asymptotic late time value) to a curved shape with decreasing temporal index a. This 
is qualitatively consistent with observations, where these type of behaviour are observed. We have demonstrated that just 
after the peak of the last pulse, the decay index of the prompt emission tail can reach very large values, far greater than the 
typical average value observed during the SDP by Swift, of 3 < a < 5 (Nousek et al. 2006). Larger values for the temporal 
index, however, are sometimes observed close to the end of the prompt emission (for example in GRB050422, GRB050803 or 
GRB050916; see figure 2 from O'Brien et al. 2006), in accord with our model. 

Because of the large number of free parameters, the fitting of actual data should be handled with care, and there may be 
various degeneracies involved. The results of such fits to data should also be taken cautiously because of the difficulty in prop- 
erly resolving distinct pulses in the prompt emission. For different reasons (such as noisy data, coarse time bins, pulse overlap, 
etc.), a group of distinct pulses may be fitted by a single broader pulse, resulting in an over-prediction of the fiux during the 
SDP, as well different spectral and temporal evolution, which might lead to a misinterpretation of the SDP. Nevertheless, when 
handled with care, a fit of our model to a good combined data set of the prompt GRB emission and its SDP tail can serve as a 
powerful test of the HLE model for the SDP, and thus help distinguish between different models for the prompt GRB emission. 

J. G. gratefully acknowledges a Royal Society Wolfson Research Merit Award. 
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Table 1. Parameters of the pulses for figure\^ (top panels) 
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Table 2. Parameters of the pulses for figure\^ (middle panels) 
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-Fpeak/ -'^O 
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0.75 


0.85 


0.95 


0.55 


0.35 
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0.11 


0.11 


1.03 



Table 3. Parameters of the pulses for figure\^ (bottom panels) 
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Figure 1. lUustration of Equal Arrival Time Surfaces (EATS). The particular EATS shown here are for a coasting shell (m = 0), and 
are ellipsoids (Rees 1966) with an ellipticity e = f3 and a semi-major to semi-minor axis ration ofT (which for display purposes is only 
3 here). The red solid lines correspond to the inner (Ro) and outer (Rf) radii of the emission region. We call Rl{T) the outermost 
radius of the EATS at observed time T . Shown are the EATS for the limiting cases corresponding to Rl = Ro (T = Tej + To; thick 
short-dashed line) and Rl = Rf (T = Toj -'rTf ; thick long-dashed line) as well as representative cases for Rl < Ro, Ro < Rl < Rf and 
Rl > Rf (dotted lines). For Rl < Ro the EATS does not intersect the emission region, and therefore the first photons start reaching 
the observer only at T = Tq + To from R = Ro along the line of sight. At Ro < Rl < Rf the flux typically rises (for AR < Ro ). At 
T = Toj + Tf the last photons from the line of .sight (at R = Rf) reach the observer, while for T > T^j + Tf the front part of the EATS, 
which would otherwise contribute the most to the observed flux, sticks outside of the emission radius resulting in a sharp decay in the 
observed flux, which is then dominated by emission from large angles relative to the line of .sight (RLE). 




Figure 2. Lightcurves of a single pulse at different normalized frequencies, vjvo. The low and high energy slopes of the spectrum are 
h\ = —0.25 and 62 = —1.25, while a = 1. AR/Ro = 1, so that Tf = 2 and Tf = 1. Left: Normalized flux density shown as a function 
ofT in linear scale. Right: Flux density shown as a function of T in logarithmic scale. 
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log(T) Tx(T„,/Tj 



Figure 3. The same pulse as in figure\^is shown for different values ofTj for {v/vq) = 0.01, 0.1, and 1 in the first, second and third 
panel respectively (in logarithmic scale). The fourth panel shows the case (I'/i'o) = I in linear scale in order to show the shape of a pulse 
having its peak before T = Tf. The normalized flux density is shown as a function of T X To^i/Tf^ where the subscript i denotes the i'th 
pulse, so that all the Tj^i would appear to coincide, and the decay times of the different pulses would appear to be the same. 
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Figure 4. Effect of the variation of a on the shape of a pulse for [v/uq) = 0.01, 0.1, and 1 in the first, second and third panel respectively 
(in logarithmic scale). We can see that increasing a makes the pulse sharper. The constant parameters are To = 1 s and Tf = 2 s and 
the low and high energy spectral slopes are bi = —0.25 and 62 = —1.25. 




T T 



Figure 5. Evolution of the temporal indexes ctgj = — dlogF^/dlogT (left panelj and Oon = —dlog Fv/dlogT (right panelj with 
normalized observed time T, at different observed photon energies (for Eg = 300 keVJ. Different line styles are used for the different 
energies. The color coding shows the temporal indexes for several values ofTf = AR/Rq: (green), 0.1 (blue), 1 (black) and 10 (red). 
The low and high energy spectral .slopes are bi = —0.25 and 62 = —1.25, while a = 1. 




Figure 6. Left panel: Comparison of the evolution of the spectral (13 - thin lines) and temporal (a^j - thick lines - and Oon - very 
thick lines) indexes at different photon energies (Eg = 300 keV). The low and high energy spectral slopes are 6i = —0.25 and = —1.25, 
while a = 1. Middle panel: Evolution of spectral index /3 over the Swift XRT energy range (umin < y < fmax with J^min = 0.2 keV 
and i^max = 10] keV. Shown are the local values of P at Vuiin ('short dashed linej, J^max/2 (dotted linej and t-max ^ong dashed linej, 
as well as the average values of jl over the XRT range, taken either over v (thin solid linej or over logi^ (thick solid linej. The low 
and high energy spectral slopes are h\ = —0.25 and 62 = —1.25, while a = 1. Right panel: Evolution of the observed spectrum with 
time. The spectrum, vF,j / (uoFq) , is shown as a function of the normalized frequency, v/vq, for different values of the normalized time, 
logjQ(r/T^), where we have used Tf = 1. The red thin lines correspond to the rising stage of the pulse (T < Ty), while the black thick 
lines are for its peak (T = Tf) and decaying stage (T > Tf). The low and high energy .spectral .slopes of the spectrum are bi = —0.25 
and 62 = —1.25, while a = 1. 
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Figure 7. For all panels, The black lines show the individual pulses, while the red line shows the total prompt emission. The normalized 
observed frequency is u/uq = 0.1 (a): Prompt emission with six pulses, all having the same following parameters: m = 0, d = —1, a = 1, 
bi = —0.25, b2 = —1.25, Tq = 2 s, Tf = AR/Rq = 3 and = 1- The ejection times Tcj are (from the first to the last pulse): 

—2 s, 8 s, 18 s, 28 s, 38 s and 48 s. (b): Same as top left panel, except for varying Fp^ak/F'o, which is from the first to the last pulse: 
0.1, 0.5, 1, 0.8, 0.3, and 0.2. (c): Same as top left panel except for varying Tq while AR/Ro = 3 remains constant, whose values are 
(from first to last pulse): 6 s, 5 s, 4 s, 3 s, 2 s, 1 s, which correspond to Tf = 24 s, 20 s, 16 s, 12 s, 8 s, 4 s. To keep tej,i = — To,i the 
ejection times in this case are: —6 s, 4 s, 14 s, 24 s, 34 s and 44 s. (d): Same as top left panel except for varying AR/Rq while keeping 
Ro and therefore To constant. The values of AR/Rq are (from first to last pulse): 5, 4, 3, 2, 1, 0. Since To = 2 s, this corresponds 
to Tf = (1 + AR/ R())To = 12 s, 10 s, 8 s, 6 s, 4 s, and 2 s, respectively, (e): Same as top left panel except for varying AR/Rq while 
keeping Rf constant and therefore Tf and Tq/Ro also remain constant, while both Ro and To vary. From first to last pulse, AR/Rq = 
10, 3, 1, 0.3, 0.1, and 0.03, and since Tf = 8 s this corresponds to To = 0.727 s, 2 s, 4 s, 6.15 s, 7.27 s, and 7.77 s. The final peak 
frequency Vp{Tf) = (To/Tf)vo at Tf is also kept constant, so that from the first to the last pulse vjvo = 0.0091, 0.025, 0.05, 0.0769, 
0.0909, 0.0971. (f): example of a more realistic prompt emission consisting of three pulses with Tej = —1 s, 13 s, 21 s, Tg = 2 s for all 
three pulses, AR/Rg = 3, 2, 1, and -Fpoak/^o = 0.7, 1, 0.7. 
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Figure 8. Comparison between the emission from several partially temporally overlapping pulses (individual pulses are in non-solid 
black lines, while the total prompt emission is in a solid red line), and a tentative fit to these pulses using a single broad pulse (thick solid 
green line). The same single broad pulse is used as a tentative fit for three different prompt emissions, with 3 (top panelsj, 6 (teiddle 
panels ) and 12 (tottom panels ) pulses. The parameters of the pulses are shown tables\^ [2| andO The normalized observed frequency is 
v/vo = 0.1 Left panels: linear scale. Right panels: Logarithmic scale. 
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APPENDIX 



APPENDIX A: DETAILED CALCULATION OF THE FLUX 



In order to calculate the flux density that reaches the observer at an observed time T, we closely follow Granot, 
Cohen- Tanugi and DoCouto e Silva 2008: we integrate over the Equal Arrival Time Surface (EATS), i.e. the locus of points 
from which photons that are emitted by the shell at a radius R, angle 6 relative to the line of sight, and a lab frame time 
reach the observer simultaneously at an observed time T. The lab frame time and the shell radius are related by 



^ dr R 

7^ ~ — + 



R 



/3c " c ' 2(m + l)cr2 
Prom simple geometrical considerations, the EATS is given by 



(1 + ^) 



R cos 9 



R 

c 



2(m+ l)r2 



Since r 3> 1 we can consider only small emission angles 9 <^1, for which cos9 1 — 9'^ /2, so that the EATS reads 



Rl 



[y'-+'+y{m + l)irL9f] 



(Al) 



(A2) 



(A3) 



(l + z) 2(m + l)cr2 

where we have introduced the normalized radius y = R/Rl, as well as Rl = Rl{T) that is the largest radius on the EATS 
at time T, and Fl = r(7?L). Since Rl is always obtained along the line of sight (at 9 = 0), 



Substituting eq. l|A4p into eq. I|A3P implies 



R^ = 2{m + l)crij^-^ 



y 



Ro 



To 



{l + z)Ro 
2(m + l)cr2 



y 



dy 



y-^ + my"-^ 



2 (m + l)r2 ' dy 2(m + l)r2 ' 
where fi = cos 9. The Doppler factor between the comoving frame and the lab frame is given by 

2r _2{m + l)rLy~ 



^^^^V " r(i-/3M) ^ i + {r9y 



i/2 



y~ 



(A4) 



(A5) 



(A6) 



Remembering the reader that T = T^j + Tf is the time at which the last photons that are emitted along the line of sight 
(from Rf and 9 — 0) reach the observer (which can be deflned here by Rl{Tc] + Tf) = Rf), from equation l|A4P its general 
value is 



In the limit AT? ^ 0, Tf ^ Tq. 



The observed flux is then obtained by integration over the EATS (Sari 1998; Granot 2005), 



F.{T) = 



(1 + ^) 
4ndl 



dLv = 



l + z 



5'^ dL',,1 — 



S-Kdl 



dy^5\y)L',,{y) , 
dy 



(AS) 



where dI/(,/ = L'^, {R)dud(f) / Atx L'^,{r)djj./2 — ^L'^,{R){djj,/dy)dy due to symmetry around the line of sight (no dependence 
of the emission on the azimuthal angle 4>), L'^,{R) is the total comoving spectral luminosity of the shell (the emitted energy 
per unit time and frequency), v' = v^l + z)/S, and dL{z) is the luminosity distance of the source. The limits of integration 
over y are 



t/min 



yn 



= min I 1 



Ro 



' Rl{T) 



1 T ^ Tej + To 



To 



min ( 1, 



Ro + AR 
Rl{T) 



-l/{m + l) 



T > T.j + To , 

T SC Tej + Tf 
T ^ Tcj + Tf 



(A9) 



For T ^ Tcj + To we have Rl{T) ^ Ro and therefore j/min = J/max = 1 and Fu{T) — 0. This is since the EATS does not 
intersect the emission region for Rl < Ro, and only touches it at one point, {R, 9) = (Ro, 0), for Rl = Ro {T = Tcj +To). The 
observed flux becomes non-zero for Rl > Ro, corresponding to T > Tcj + To. Substituting eqs. I|A5P and l|A6ll into eq. I|A8P 
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finally gives 

Al Power-law spectrum 

While a single power law emission spectrum is not very realistic, it already shows many important properties that also 
appear for a Band function emission spectrum (considered in the main text). This is the reason why this case is described 
here. The luminosity is then 

^ (1)° iQ' ^ (srfF^) ' ■ 

where the comoving spectral luminosity also scales as a power law with radius when the emission is over a finite range of 
radii, Vq is a fixed frequency in the comoving frame. 



A 1.1 Emission from an infinitely thin shell at radius iio 



We first study the case where the whole emission comes from a single radius -Ro, 



L'„, 



= io 



Ro5{R - Ro) = L'o 
m + 1 



y 



mb/2 



yi-ninS{y - yn 



where this is valid only for T ^ 1 that corresponds to Rl ^ Ro, for which y, 
then implies 

b 



ymin5{y - ymin) 

(A12) 

, = Ro/Rl = [(r-rcj)/ro]-'/''"+'>. Eq. 



F.{T ^ Tej + To) = ^i^Lo 



47rdi 



There are two times of particular relevance here: the radial time Tr{Ro) = To = (1 + z)Ro/[2c{m + l)ro], which is the time 
past Tej when the first photons start reaching the observer, and the angular time Te{Ro) = (1 + 2:)-Ro/(2cro) — {m + l)Tr{Ro) 
that sets the time-scale for the width of the pulse. One can rewrite the expression for the observed flux density as 



F^{T > Tej + To) 



(1 + ^) 
4ndl 



Lo 





'T-Ts' 




Te{Ro)_ 



(A14) 



where Ts = Tcj + To — Te{Ro) = Tcj — rriTo is the reference time for the power-law flux decay of the pulse, and is exactly 
Tg{Ro) before the onset of the pulse. Since the emission itself occurs at one particular radius (-Ro) it depends only on the 
Lorentz factor at that radius radius, and is independent of m. In particular, re(-Ro) = {m + l)To and the pulse peak flux are 
independent of m. The value of m affects only the onset time of the pulse (T — Tcj + To) and the reference time for the power 
law flux decay. For internal shocks we expect a coasting shell (m = 0) for which Ts — Tcj and Tb{Ro) = To. It can easily be 
seen that the HLE relation, a = 2 + P where Fi, oc T~"v~^ , is satisfled here as /3 = — 6 and a = 2 — h— 2 + 13. 



A 1.2 Emission from a region of finite width 

We now turn to the case where the emission comes from a range of radii between Ro and Rf = Ro + AT? > _Ro- The 
comoving spectral luminosity in this case is L'^i — L'o 

{R/RoT W /v'o) , and the flux density is given by (Granot, Cohen- Tanugi 

& do Couto e Silva 2008): 

which, for internal shocks (m — 0) becomes: 

T.(T > Tej +To) = (1±|) l^ii^T-^ [min(T,f,)^+-' - l] . (A16) 

It is therefore obvious that for T ^ Tcj + Tf the HLE relation is valid, where the reference time is the ejection time Tcj , as 
in this case the spectral slope \s (5 — —h and the temporal slope is a = 2 — 6 = 2 -I- /3. In this sense a flnite range of emission 
radii with m = is similar to emission from a single radius, as in both cases the HLE relation a = 2 -|- /? is strictly valid 
immediately from T ^ Tej -|- T/, for some reference time, though in the latter case the reference time for which this is valid is 
equal to the observed ejection time only for m = 0. For emission from a flnite range of radii with m 7^ the relation a = 2 + /3 
is approached asymptotically at T — Tcj 2> T/ . 
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A2 Band function spectrum: general case and late time dependence 

In the main text we have given the flux in the speciflc case of internal shocks, m = and d = — 1. We derive here the 
flux for any values of the parameters m and d. 



A3 Emission from a single radius 

When the whole emission comes from a single radius Ro, the peak frequency is v'^ = v'q, and the luminosity is thus 



L'^i = L'f^S yj^j Ro5(R - -Ro) = L'qS ^JJi j yminS{y — J/min) 

Using this luminosity (eq. |A17) 1 in the integral for the flux (eq. |A10| 1 results in 



F^{T > Toj + To) = ^i^Lo 



T-T, 



Te{Ro) 



S 



fo 



T-T, 



TeiRo) 



(A17) 



(A18) 



where as in § lAll Ts = Tej — mTo is the reference time for the power-law flux decay of the pulse, Te{Ro) — (m + l)ro = 
(1 + z)Ro/2crQ is the angular time at Ro, and Ep{T) = hVp{T) is the photon energy corresponding to the peak of the Band 
function spectrum. Note that for m = 0, = Tcj. One can express the argument of S as 



E 



i^p{T) Ep{T) vo 

Reminding that Fs 
observed flux as: 



T-Ts 



Tg{Ro) 



fo 



E 



(A19) 



La(l + z)/{'^-Kd\), we then use the explicit expression for the Band function (eq. [T]) to express the 



F^jT > Tcj +ro) 



f [ 



= < 



rfl(«o)J 

r T-T, 1' 



"a ) 



Ta(Ro) 



(A20) 



A4 Emission from a range of radii 

In the case of emission with a Band function spectrum over a flnite range of radii, Ro < R < Rf = Ro + AR, we remind 
that the comoving luminosity is: 

L'.'-L'J^ys(-;^) (A21) 



^RoJ \i^UR) 

Introducing l|A2ip into l|A10P we obtain the general expression of the flux 



Ft, -a + m/2 P"^^ a-l-m/2 / m+1 \^ f 1^ 



ir = yZr'-J '^^^'^ '"^ U+V--^ ) ^(-) (A22) 



where y = y/y^in, y^ax = min[y^;i^„, (Tj/ro)^''(™+i'], and 

/ / \ m/2 — d / —fri — l\ /™ i —m — 1 m— 1 \ 

ii = ii f ^) ( ^+y^^ ) = ii r /^-^ ( — ) , (A23) 

t'p V ymin / \ m + l J UO \ m+1 ) 

and the expression for v' jv^ assumes that v'^ — v'o{R/Ri)Y = J^o{y / ymiuY ■ At late times, T — Tcj S> T/, we have y <^ 1, 
j/max = Tf/To, and v'/i^'p ~ (i^ /'^o)y~^~''~"'^^y^]^^^^ / {m + 1) increases with time so that S{v' /Up) oc [v /v'p)''^ and F^ oc 
{v/uof^y[^^^^^^'^~^'''' = {u/voP[{T - ^cj)/^o]''2~^ i.e. the HLE relation q = 2 + /? is satisfled. 



In the case for internal shocks, with m = and d — —1, v' /v'p becomes independent of y and can be taken outside the 
integral (y' jv'p — {v / vo) / y-mi-n — {u/vo)T), leading to the much simpler expression of the flux seen in the main text (eq. [T}. 
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APPENDIX B: EVOLUTION OF THE TEMPORAL AND SPECTRAL INDEXES 



This appendix explicits the evolution of the temporal and spectral indexes with time. 
BO.l Single emission radius 

Where the luminosity is a delta function with radius at radius -Ro, we obtain 

-bi+f{l + bi)u/uo f^XbUo/u 



OLoTl 



-&2 T ^ XbV(,/v 

(2 - 6i)f /(I + T) + f (1 + W)v/vo T < Xbvo/v - 1 
(2-fe2)T/(l + T) f';tXbVo/u-l 



(Bl) 



(B2) 



(B3) 



2-&2 



T ^ XiVo/v 



We then have a very simple relation between a^j and (3: a^j = 2 + f3, as expected at asymptotically late times for HLE, 
just that for aej it is satisfied all along for the local values of the temporal and spectral indexes. At late times Oon approaches 
Qcj and a similar relation approximately holds between aon and /3 {cton ~ 2 + f3). 

BO. 2 Emission from a finite range of radii: Ro < R < Rj 

In this case, the spectral index /3 is still given by eq. IjBip . while the two temporal indexes are: 

' (2 - &i) + r(l + bi)^ - (2 + a) ffT)fil-°i f <mm{AR/Ro,xt,yo/u-l) , 



aon = < 



i2-h)jj^ + Til + b,)^ 



{2-b,)j^-{2 + a) ^^['+ll:^ 



which limits at very early and very late times are 

-1 r<i, 



AR/Ro <T < xtuo/u - 1 , 
Xbi^o/u - 1< f < AR/Ro , 
T > ma,x{AR/Ro,xtuo/i' - 1) 



(B4) 



(B5) 



and 



2 - &2 T > 1 , 

2 - bl - (2 + a)/(l - f ^ ^ bi)iy/,^Q f < mm{Rf/Ro, xtuo/ry) ., 

2-6i+T(l + 6i)i//t'o Rf/Ro<f<XiUo/u, 

2 - 62 - (2 + a)/(l - f -"-2) XbVo/v <f < Rf/Ro , 

2 — 62 T > ma.x{Rf/Ro,Xbiyo/i^) 

According to equations ijBip and ljB6P Qej has a simple relation with /3: 



(B6) 



/3 + 2- (2 + a)/(l -f'-"-^) f< Rf/Ro {T < T 



f) . 



(B7) 



/3 + 2 



f> Rf/Ro {f>ff) 



Note that in the limit T ^ (T — > 1), at very early times, just after the onset of the spike, tJej —00 while Ofon — 1. 
Moreover, the simple HLE relation, aej = 2 + /3, is valid as soon as T > T/, for any value of Tf. This is a relation between 
the local values of aej and /3, that hold as both change with time, and is strictly valid from T > Tf only for m = and 
d = —1. For general values of m or d this local HLE relation would be valid only at late times, T S> Tf - Note, however, that 
for alternative other definitions of the temporal index, such as aon, this relation is only approached at late time: aon ~ 2 + (3 
for f > Tf and f » 1. 
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APPENDIX C: EXPONENTIAL TURN-OFF OF THE EMISSION WITH RADIUS 



Throughout the paper we have assumed that the emission abruptly turns off at i?/. This results in a sharp change in the 
temporal index at T/, which usually corresponds to a sharp peak for the pulses in the prompt GRB light curve. Observations 
sometimes show pulses with a round peak, which may be hard to fit with spiky theoretical spikes. Such rounder peaks for the 
pulses may be obtain within the framework of our model by introducing a more gradual turn-off of the emission at R > Rf. 
For convenience, we parameterize this here by assuming that the luminosity starts decreasing exponentially with radius at 
R > Rf. For simplicity we consider here only AR > 0, but the results are similar for AR = 0. Similarly, only the case for 
internal shock (m = 0, d = — 1) is considered here. Thus, we introduce the following comoving spectral luminosity: 



^°(^)"^(^) Ro^R^Rf 

(CI) 

^'4^ysU)e'^ R>Rf 



where q is the decay constant (a larger q corresponds to a sharper turn-off of the emission). 

For I ^ T ^ Tf the observed flux is identical to that without introducing the gradual emission turn-off, and is therefore 
given by eq. 



F,(f ^ ff) = Fof- 



- - \2-l-a 

<T,Tf)] 



S[—T] . (C2) 



The flux for T > T/ is obtained by calculations very similar to those of section 13. 2| and reads 

F^f^ff) = Fof~''S{fu/po)[ff+'' -l + .J{f)'j , (C3) 



J(T) = (2 + a) / dyy" e"^""^^'^^ , (C4) 

were Q = AR/(qRo), and we remind the reader that y = Ty for m = 0. The expression for the flux is thus very similar its form 
for an abrupt turn-off of the emission at Rf, but with the additional term J{T) that adds some flux at T > T/ (representing 
the added contributions from R > Rf). For a = 1 we have 

J{f,a = 1) = 6Q^ + 6Q^ff + 3Qff - e~'-^~^f^^'^ (gQ'-' + 6Q'^f + 3Qf^^ . (C5) 
At late times J{T, a = 1) approaches a constant value, 

J^^6Q' + 6Q'ff+3Qf]r.l ,2 (C6) 

where we have replaced Q and Tf by their dependence on q and AR/ Rq. Since J{T) appears in eq. I|C3P in a sum with f^+^-l, 
it will dominate the observed flux at late times if Joe > T^^" — 1 or equivalently if g < gcrit where Joo(gcrit) = Tj^" — 1. 

The left panel of flgure [CTI shows qcrit as a function of Tf — AR/Rq for a = 1, and it can be seen that the limiting values 
of (jcrit are 1 for ff < 1, and (7 + 22^''^Y'^ + 1 + 3/(7 + 22^''^ « 4.59 for fj > 1, so that (?crit is always of order unity. 
Therefore, for q <C gcrit ~ 1 the late time flux is dominated by contributions from R > Rf, the peak of the pulse is rounder 
and the peak flux is higher compared to an abrupt turn-off of the emission with radius, which is approached in the opposite 
limit of g S> gcrit ~ 1. This can nicely be seen in the right panel of flg. IClj which shows the shape of a pulse for different 
values of g, including the limiting case of g ^ cx), which corresponds to an abrupt turn-off of the emission at Rf. 

Such an exponential turn-off could therefore be useful when fltting our our model with data, in order to reproduce 
round-peaks pulses. Of course, one should be aware that this adds a free parameter (g or Q), and might thus increases the 
degeneracy between the different flt parameters. Therefore, adding this extra model parameter should be done only when it 
is required by the data. 
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Figure CI. Left: Dependence of the critical decay index qcrit of the exponential cut-off on Tf = AR/Rq, in semi-logarithmic scale. It 
is defined by Joo(gcrit) = — 1, i-e. at late time, the contribution to the flux from the exponential cut-off is equal to the contribution 
from the emitting region between Rq and Rf . Right: Comparison of the shape of pulses with and without the exponential turn-off of the 

luminosity for a ratio AR/Rq = 1 in logarithmic scale. The solid line shows the shape of the pulse for an abruptly stopping luminosity 
(no exponential turn-off), the other lines show the pulse shape for different values of the decay constant q = 0.1 of the exponential 
turn-off. 



